Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1204 1489 (N25W07 L=319)   478    119        0 (1178/1430-1224/1550) -60.3 W04-W08 N21-N28
      2  460 1389 (N12E41 L=272)  1069    189        3 ( 416/1328- 510/1464) -41.6 E45-E36 N08-N16
      3  484 1264 (N04E37 L=276)   515    164        1 ( 426/1190- 547/1332)  53.1 E39-E32 N00-N08
      4 1671 1184 (N12W40 L=353)  1978    416        6 (1554/1068-1817/1264) -10.3 W31-W53 N06-N16
      5 1146 1189 (N05W05 L=318)   390    120        4 (1105/1153-1194/1214)  38.7 W03-W08 N03-N07
      6  624  981 (S11E24 L=288)   278     87        1 ( 577/ 964- 664/1003)  18.4 E27-E21 S12-S09
      7  155  827 (S23E64 L=249)   809    159        4 ( 125/ 797- 185/ 856) -31.0 E67-E60 S24-S20
      8 1568  827 (S11W37 L=349)   368     96        0 (1532/ 800-1613/ 851)  -8.1 W33-W40 S11-S09
      9 1198  669 (S25W15 L=328)  1701    389        1 (1016/ 576-1389/ 796)  39.0 W04-W25 S32-S14
     10  876  540 (S37E05 L=308)   609    140        2 ( 814/ 514- 951/ 567)  24.0 E09-W00 S39-S34
     11 1571  404 (S36W52 L=  4)   333     85        0 (1540/ 385-1597/ 429)  -8.1 W47-W54 S36-S34
  

    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