Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  877 1507 (N40W03 L=123)   500    105        1 ( 842/1463- 916/1553)  29.7 W00-W07 N36-N42
      2 1662 1510 (N17W58 L=177)   452    105        2 (1641/1480-1693/1529)   0.0 W54-W61 N15-N18
      3  303 1367 (N43E48 L= 71)   525     78        1 ( 279/1327- 331/1395)  63.4 E51-E45 N40-N45
      4  178 1167 (N31E60 L= 59)   872    130        1 ( 154/1124- 207/1200) -59.0 E62-E57 N28-N33
      5 1048 1070 (N09W02 L=122)   473    107        1 (1031/1010-1065/1137)  76.0 W00-W03 N05-N13
      6  624 1044 (N17E24 L= 96)   244     46        0 ( 598/1028- 658/1062) -36.9 E25-E21 N15-N18
      7  202  957 (N18E60 L= 59)   565     75        0 ( 184/ 912- 219/ 996) -78.7 E61-E58 N15-N20
      8 1046  372 (S34E16 L=103)   641     90        1 ( 997/ 350-1093/ 399) -39.8 E19-E13 S36-S31
  

    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