Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1168 1762 (N44E18 L=196)   614    156        0 (1107/1699-1215/1822) -21.8 E26-E11 N40-N47
      2  324 1446 (N00E60 L=154)  1016     87        3 ( 299/1399- 355/1482) -76.0 E64-E56 S02-N02
      3 1178 1377 (N18E02 L=213)   773    182        0 (1107/1300-1240/1430) -21.8 E06-W04 N14-N19
      4  360 1341 (S05E51 L=163)   473     64        2 ( 335/1308- 383/1378)  90.0 E53-E49 S07-S02
      5  211 1070 (S24E60 L=154)   531     52        1 ( 192/1039- 229/1105)  90.0 E62-E57 S25-S21
      6  768  910 (S20E11 L=203)  1193    433        0 ( 579/ 816- 947/ 969)  41.8 E21-E02 S29-S11
      7  396  808 (S36E37 L=177)   747    180        1 ( 363/ 731- 455/ 870)  61.4 E40-E32 S41-S30
      8  893  754 (S25W00 L=214)   590    177        1 ( 826/ 711- 982/ 783)  56.3 E02-W04 S29-S20
      9  964  628 (S30W09 L=223)   956    284        5 ( 823/ 610-1096/ 656)  29.4 W00-W16 S34-S25
     10 1256  502 (S28W32 L=246)   330     84        0 (1216/ 482-1297/ 526)   9.5 W28-W34 S28-S26
     11 1559  459 (S20W56 L=270)  1007    149        0 (1489/ 420-1623/ 503)  56.3 W52-W59 S24-S15
     12  379  458 (S62E47 L=168)  1057    149        2 ( 317/ 429- 434/ 490)  -2.6 E60-E35 S64-S59
  

    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