Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1232 1377 (N21W22 L=202)   216    131        0 (1181/1336-1279/1421)  12.5 W17-W26 N20-N22
      2 1400 1363 (N16W32 L=211)   395    132        3 (1345/1320-1442/1403) -74.1 W29-W34 N12-N19
      3  979 1358 (N27W06 L=186)   392    160        0 ( 926/1310-1051/1384)   0.0 W01-W10 N25-N28
      4  480 1297 (N36E31 L=149)  1373    504        1 ( 337/1205- 675/1355)  -3.9 E44-E15 N31-N38
      5 1796 1268 (S01W59 L=238)   560     68        2 (1763/1235-1829/1287) -15.9 W55-W62 S03-N00
      6 1431 1185 (N05W28 L=207)   420    103        0 (1397/1139-1466/1240)  29.7 W24-W31 N03-N07
      7 1460 1019 (S04W25 L=204)   342     83        2 (1439/ 974-1481/1052) -90.0 W23-W26 S07-S02
      8 1864  955 (S21W58 L=237)  1154    281        1 (1798/ 825-1929/1040)  21.8 W48-W69 S27-S17
      9 1354  664 (S22W10 L=190)   757    293        0 (1236/ 593-1429/ 763)  17.4 W01-W17 S24-S18
     10  502  472 (S14E50 L=130)   310     54        0 ( 481/ 445- 524/ 500) -90.0 E50-E48 S15-S11
  

    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