Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  426 1394 (N30E45 L=194)  1411    313        2 ( 346/1318- 551/1450)  16.4 E51-E32 N25-N34
      2  992 1398 (N28W00 L=239)   371    173        1 ( 905/1365-1050/1420) -18.4 E06-W03 N26-N30
      3  263 1377 (N29E63 L=176)   645     83        1 ( 225/1352- 304/1394)  10.3 E68-E57 N27-N30
      4  476 1321 (N26E38 L=201)   673    145        1 ( 410/1301- 540/1346)  10.3 E43-E32 N24-N27
      5  203 1184 (N17E65 L=174)  1775    240        4 ( 179/1086- 229/1290) -69.0 E69-E60 N10-N23
      6  848 1220 (N18E10 L=229)   222     99        2 ( 825/1193- 887/1269)  45.0 E11-E07 N16-N20
      7 1316  935 (S03W18 L=257)   669    172        3 (1243/ 882-1358/ 997)  35.5 W12-W20 S05-N01
      8 1803  914 (S08W56 L=295)   872    121        2 (1780/ 857-1828/ 966)  56.3 W53-W59 S11-S05
      9 1743  857 (S11W50 L=289)   505     62        1 (1724/ 825-1767/ 889)  76.0 W48-W52 S13-S09
     10  774  565 (S23E20 L=219)   776    202        2 ( 752/ 482- 806/ 651)  80.5 E21-E16 S29-S17
     11  539  550 (S24E38 L=201)   357    104        0 ( 491/ 525- 581/ 582)  18.4 E42-E33 S25-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