Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1608 1565 (N45W26 L=142)   375     55        0 (1578/1543-1633/1591)  -9.5 W22-W28 N42-N46
      2  925 1465 (N15E18 L= 98)   494    154        1 ( 899/1393- 957/1539) -45.0 E21-E13 N11-N19
      3 1052 1434 (N17E09 L=106)   302     97        0 (1013/1420-1101/1447)  21.8 E11-E06 N15-N18
      4 1144 1354 (N14E02 L=114)   307    112        0 (1126/1295-1164/1403) -51.3 E04-W00 N11-N16
      5 1828 1090 (N21W52 L=168)   749    157        0 (1814/1012-1844/1163) -63.4 W49-W55 N16-N25
      6 1343 1039 (N02W17 L=133)   261     33        1 (1314/1021-1369/1060)  14.0 W15-W19 N00-N03
      7 1594  880 (N01W37 L=153)   444     52        0 (1562/ 853-1620/ 909)   0.0 W35-W40 N00-N02
      8  862  585 (S36W03 L=119)  1153    277        1 ( 742/ 552- 999/ 630)  33.7 E05-W10 S41-S29
      9 1062  476 (S37W21 L=136)   444     87        1 (1043/ 435-1083/ 530) -56.3 W18-W23 S39-S33
     10 1488  433 (S25W51 L=167)   734    127        1 (1443/ 382-1535/ 493) -77.5 W49-W53 S29-S20
  

    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