Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  611 1636 (N30E37 L=192)   362     60        0 ( 579/1618- 643/1652) -11.3 E40-E34 N28-N30
      2  817 1492 (N21E18 L=211)   637    235        1 ( 724/1447- 932/1555) -17.4 E26-E10 N19-N24
      3  445 1249 (N03E40 L=190)   409     64        0 ( 424/1206- 468/1291)  90.0 E40-E38 N00-N05
      4  717 1128 (S02E19 L=210)   424    137        1 ( 658/1091- 791/1167) -12.5 E23-E14 S03-N00
      5  796  918 (S14E13 L=216)   200     44        1 ( 765/ 905- 821/ 936)   0.0 E14-E11 S14-S12
      6 1377  809 (S14W24 L=253)  1340    294        1 (1253/ 721-1517/ 934)  49.1 W17-W31 S20-S05
      7 1111  829 (S16W07 L=237)   207     10        0 (1092/ 806-1129/ 848)  56.3 W06-W08 S17-S14
      8  626  780 (S24E23 L=206)   253     64        0 ( 595/ 765- 665/ 793)  11.3 E25-E20 S24-S22
      9  779  690 (S28E13 L=216)   464    147        0 ( 701/ 643- 838/ 733)  31.0 E18-E08 S31-S24
     10 1564  654 (S21W41 L=270)   260     84        0 (1535/ 610-1591/ 698)  80.5 W39-W41 S23-S17
     11 1031  364 (S48W11 L=240)   272     54        0 ( 998/ 350-1069/ 379)  21.8 W07-W13 S49-S46
     12  645  348 (S55E27 L=202)   344     42        1 ( 621/ 330- 667/ 365)  21.8 E29-E24 S56-S52
  

    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