Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  820 1686 (N30E34 L=231)  1259    250        2 ( 742/1615- 896/1741) -17.4 E41-E25 N26-N33
      2  645 1647 (N23E43 L=221)  1378    210        3 ( 567/1558- 714/1725)  85.9 E45-E40 N16-N30
      3 1323 1408 (N26W08 L=272)  2216    490        4 (1209/1194-1382/1590) -49.4 E04-W15 N13-N35
      4 1225 1194 (N10W07 L=271)   213     75        1 (1196/1163-1254/1225) -21.8 W04-W09 N09-N11
      5 1634 1036 (N13W37 L=301)   415    100        0 (1618/ 991-1647/1099) -59.0 W34-W38 N09-N15
      6  291 1054 (S19E50 L=215)   634    124        1 ( 244/1008- 341/1082)  36.9 E53-E44 S22-S16
      7  849  996 (S10E09 L=255)   201    123        0 ( 827/ 944- 888/1028)  90.0 E10-E07 S13-S07
      8  659  912 (S19E19 L=245)   467    211        0 ( 625/ 856- 706/ 974) -39.8 E21-E14 S21-S14
      9  288  914 (S27E48 L=216)   462     92        0 ( 243/ 904- 341/ 931)  10.3 E53-E42 S28-S26
     10 1744  766 (N01W53 L=318)   317     36        0 (1722/ 746-1770/ 782)  33.7 W51-W55 N00-N02
     11  616  451 (S50E14 L=251)  1701    266        3 ( 511/ 380- 724/ 538) -12.3 E24-E01 S54-S46
     12 1220  485 (S29W29 L=293)   254     67        0 (1193/ 460-1255/ 509)  63.4 W27-W30 S31-S26
     13 1125  190 (S52W48 L=312)   388     42        0 (1110/ 168-1142/ 217) -16.7 W43-W53 S53-S50
  

    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