Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1037 1611 (N44W12 L=345)   277     49        1 (1010/1585-1063/1633)  36.9 W09-W14 N41-N45
      2 1442 1245 (N13W30 L=  3)   610    173        0 (1414/1167-1486/1317) -84.3 W28-W31 N07-N17
      3  302 1169 (N23E51 L=282)   347     71        1 ( 281/1137- 324/1206) -51.3 E53-E48 N20-N25
      4  673 1139 (N17E21 L=312)   309    103        0 ( 639/1096- 708/1180) -56.3 E22-E18 N14-N20
      5  732 1019 (N09E18 L=315)   301    101        1 ( 697/ 994- 775/1046) -21.8 E20-E15 N07-N11
      6 1684  906 (S12W42 L= 15)   390    113        0 (1639/ 875-1707/ 949)  26.6 W38-W44 S13-S09
      7 1631  729 (S22W37 L= 10)   497    142        1 (1591/ 667-1664/ 781)  45.0 W34-W40 S25-S19
      8 1441  691 (S21W22 L=355)   343    126        0 (1415/ 633-1462/ 748)  60.3 W20-W24 S24-S17
      9 1108  481 (S30E03 L=330)   365    107        0 (1087/ 430-1143/ 519)  56.3 E04-W00 S33-S26
     10  414  416 (S27E61 L=272)  1656    164        4 ( 390/ 336- 438/ 496)  29.5 E75-E52 S33-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