Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  496 1565 (N53E35 L= 73)   576     82        1 ( 479/1519- 518/1606) -54.5 E38-E32 N49-N56
      2  740 1573 (N48E08 L= 99)   810    126        3 ( 684/1541- 803/1601) -31.0 E13-E03 N44-N51
      3  714 1296 (N29E15 L= 93)   366    134        0 ( 676/1246- 769/1340)  29.7 E17-E10 N26-N31
      4  219 1136 (N27E60 L= 47)   995    158        2 ( 190/1065- 256/1202)  48.8 E63-E56 N22-N30
      5  345 1126 (N25E46 L= 62)   613    171        2 ( 295/1058- 386/1186) -48.4 E51-E42 N19-N28
      6  151 1053 (N22E70 L= 37)   888    124        0 ( 129/ 999- 176/1115)  49.4 E74-E66 N18-N25
      7  819  836 (S00E16 L= 92)   723    209        2 ( 755/ 746- 876/ 909) -70.0 E18-E13 S06-N05
      8 1099  754 (S10E01 L=107)   247     69        1 (1071/ 719-1129/ 780)  18.4 E03-W00 S11-S08
      9 1734  725 (S28W44 L=152)   367     36        0 (1708/ 706-1752/ 749)  14.0 W41-W46 S29-S26
     10  481  382 (S24E58 L= 50)  2019    269        1 ( 386/ 294- 578/ 441) -58.4 E65-E52 S32-S18
  

    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