Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1375 1555 (N33W17 L=214)   284     92        0 (1338/1529-1398/1597)  56.3 W13-W18 N30-N36
      2  333 1523 (N18E59 L=138)  2482    156        5 ( 264/1461- 393/1589)  31.0 E65-E53 N13-N22
      3 1094 1513 (N25E02 L=195)   752    170        2 (1057/1446-1124/1571) -58.0 E05-W00 N20-N28
      4  351 1426 (N12E52 L=145)   751     76        1 ( 316/1390- 379/1455) -33.7 E55-E49 N09-N13
      5 1170 1356 (N16W05 L=202)   405     97        2 (1121/1324-1220/1388)   0.0 W01-W07 N14-N17
      6 1779 1278 (N21W50 L=248)  2206    342        3 (1656/1144-1864/1387) -28.8 W39-W59 N14-N26
      7  524 1233 (N01E34 L=164)   283    111        0 ( 494/1208- 573/1272) -18.4 E36-E30 N00-N03
      8  938 1198 (N03E07 L=190)   371    117        0 ( 902/1156- 977/1239) -31.0 E09-E04 N01-N05
      9  397 1121 (S07E42 L=156)   287     88        0 ( 372/1084- 436/1146)  38.7 E43-E38 S08-S04
     10 1379  902 (S08W23 L=221)   252     58        0 (1356/ 875-1409/ 927) -36.9 W21-W25 S09-S06
     11  676  802 (S23E19 L=178)   474    119        1 ( 619/ 782- 739/ 827)   0.0 E23-E15 S24-S22
     12 1310  745 (S18W22 L=219)   685    271        1 (1188/ 697-1433/ 787)  28.1 W14-W29 S22-S14
     13 1666  402 (S32W69 L=267)  2426    302        1 (1541/ 346-1782/ 488)  39.8 W58-W78 S38-S23
  

    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