Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1587 1669 (N24W65 L=180)  1315    189        3 (1492/1627-1648/1697) -21.8 W56-W71 N19-N27
      2 1459 1553 (N22W45 L=160)  1602    394        1 (1306/1513-1646/1606) -36.9 W36-W57 N14-N29
      3 1295 1336 (N15W24 L=140)   308    119        0 (1240/1322-1362/1344) -23.2 W21-W28 N13-N17
      4 1587 1286 (N04W41 L=156)   382    107        0 (1556/1235-1617/1335) -81.9 W40-W41 N00-N07
      5  983 1285 (N21W05 L=120)   300     89        0 ( 940/1269-1029/1298) -31.0 W02-W07 N19-N22
      6  428 1163 (N28E35 L= 80)  1129    249        3 ( 327/1112- 536/1217) -26.6 E44-E26 N24-N33
      7  753 1075 (N15E14 L=101)   921    371        3 ( 611/1037- 919/1105) -17.5 E22-E03 N11-N19
      8  250 1022 (N24E54 L= 62)   413     95        1 ( 227/ 986- 272/1064)  59.0 E56-E51 N21-N26
      9 1465  692 (S26W17 L=133)  1475    366        2 (1388/ 530-1500/ 860)  54.2 W09-W23 S34-S16
     10  529  368 (S21E55 L= 60)  1838    340        3 ( 401/ 340- 690/ 399) -12.5 E70-E42 S24-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