Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  946 1795 (N52E25 L=222)   580    103        1 ( 898/1779- 998/1814)   5.7 E29-E19 N50-N53
      2  314 1509 (N20E63 L=184)  1394    222        3 ( 281/1463- 357/1568) -45.0 E67-E56 N17-N23
      3  410 1381 (N13E48 L=199)   330    146        0 ( 370/1340- 447/1407) -18.4 E52-E43 N11-N14
      4  194 1218 (N01E66 L=181)   511     57        0 ( 179/1190- 209/1242) -11.3 E68-E63 N00-N02
      5 1244 1181 (N13W11 L=258)   431    207        1 (1150/1166-1334/1218)   0.0 W05-W17 N11-N14
      6 1348  951 (N00W21 L=267)  1133    238        2 (1312/ 864-1407/1063) -70.0 W18-W25 S04-N07
      7  312  862 (S19E48 L=199)   490    145        0 ( 276/ 822- 368/ 912) -32.0 E51-E42 S21-S16
      8  404  869 (S18E40 L=207)   370     78        3 ( 379/ 832- 434/ 900)  78.7 E41-E37 S20-S15
      9  814  877 (S11E11 L=236)   432    128        0 ( 746/ 861- 885/ 889)  12.5 E15-E06 S13-S09
     10 1054  795 (S13W05 L=252)   730    287        1 ( 947/ 714-1179/ 875)  45.0 W00-W12 S19-S07
     11 1430  680 (S15W31 L=278)  1244    159        4 (1354/ 569-1494/ 782)  75.1 W28-W33 S23-S08
     12 1228  497 (S30W23 L=269)   716    152        3 (1168/ 458-1317/ 532)  41.2 W18-W28 S33-S26
  

    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