Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  741 1733 (N33E42 L= 93)   487     52        0 ( 716/1708- 765/1765) -78.7 E43-E39 N31-N36
      2 1182 1686 (N43E09 L=126)  1326    217        3 (1139/1601-1239/1765) -35.5 E15-E01 N37-N48
      3  771 1455 (N16E26 L=109)  1037    342        2 ( 731/1291- 840/1591) -76.0 E29-E22 N05-N25
      4  488 1479 (N11E46 L= 89)   956    183        2 ( 452/1392- 541/1540)  84.3 E49-E44 N05-N15
      5  169 1240 (S09E70 L= 65)  2650    414        8 ( 120/1149- 206/1308) -73.3 E78-E63 S14-S04
      6  666 1169 (S03E24 L=111)   303     96        2 ( 641/1135- 684/1217) -45.0 E26-E22 S05-N00
      7 1186 1162 (N09W06 L=141)   207     53        1 (1159/1143-1214/1189) -14.0 W03-W07 N07-N09
      8  447  850 (S27E33 L=102)   301    111        1 ( 406/ 831- 472/ 893)  59.0 E36-E31 S28-S23
      9  935  861 (S14E01 L=134)   186     42        2 ( 908/ 845- 958/ 879)  45.0 E02-W00 S15-S12
     10  617  654 (S35E18 L=117)   278     55        2 ( 593/ 625- 640/ 682) -56.3 E19-E15 S36-S33
     11 1552  587 (S15W46 L=181)   485    173        3 (1488/ 545-1586/ 617)  20.6 W40-W49 S16-S12
     12 1281  325 (S38W41 L=176)  1853    454        2 (1093/ 296-1479/ 369)  26.6 W27-W53 S42-S29
  

    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