Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  485 1618 (N19E55 L= 17)   439    103        0 ( 455/1584- 519/1659)  80.5 E56-E53 N16-N22
      2 1299 1538 (N35W04 L= 77)  1315    256        2 (1238/1428-1346/1632) -47.5 E02-W09 N28-N40
      3  384 1489 (N09E56 L= 16)   831    135        1 ( 359/1442- 409/1551) -54.5 E58-E52 N06-N13
      4 1609 1165 (N19W33 L=106)  1117    257        2 (1564/1054-1668/1255)  81.9 W30-W37 N12-N26
      5 1830 1038 (N18W56 L=129)  1420    178        0 (1804/ 961-1860/1134) -56.3 W52-W60 N14-N23
      6 1314 1036 (N04W16 L= 89)   345     82        0 (1275/1010-1347/1062)  53.1 W14-W17 N02-N06
      7 1905  866 (N11W72 L=144)   607     70        0 (1875/ 849-1926/ 898)  38.7 W66-W76 N09-N13
      8  986  743 (S20W05 L= 77)   267    121        0 ( 941/ 720-1044/ 758)  45.0 W02-W07 S22-S17
      9  543  707 (S34E24 L= 48)  1017    221        2 ( 456/ 675- 639/ 747)  15.9 E31-E17 S35-S29
     10 1595  558 (S15W51 L=124)   449    101        1 (1561/ 519-1636/ 595)  51.3 W49-W53 S17-S12
     11 1558  404 (S25W59 L=132)   608    111        3 (1520/ 370-1598/ 447)  80.5 W56-W63 S27-S21
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y