Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  725 1807 (N44E37 L= 26)   709     97        2 ( 696/1778- 759/1838)  45.0 E41-E34 N40-N46
      2  884 1511 (N23E15 L= 49)   228     17        0 ( 863/1494- 905/1531)  63.4 E15-E13 N21-N24
      3 1462 1291 (N16W24 L= 88)   263     49        0 (1431/1263-1487/1314) -33.7 W22-W26 N14-N16
      4  233 1106 (S07E56 L=  8)   602    143        2 ( 192/1064- 266/1159) -56.3 E60-E52 S09-S03
      5  462 1028 (S10E36 L= 28)   337     86        0 ( 440/ 987- 479/1077)  78.7 E37-E34 S12-S07
      6  223 1029 (S12E56 L=  7)   425     76        2 ( 205/ 993- 244/1064) -63.4 E58-E54 S13-S09
      7 1651  848 (S08W43 L=107)  1319    271        1 (1564/ 809-1739/ 899)   0.0 W35-W50 S09-S04
      8  386  782 (S26E42 L= 22)   559    174        0 ( 358/ 727- 440/ 851) -50.2 E43-E37 S29-S21
      9 1087  743 (S22W07 L= 71)   215     67        0 (1063/ 704-1107/ 776)  78.7 W06-W08 S24-S19
     10 1475  579 (S26W37 L=101)   417     67        0 (1440/ 549-1503/ 618)  68.2 W34-W38 S28-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