Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1714 1474 (N20W62 L=223)  2534    379        5 (1623/1411-1809/1533) -24.0 W57-W70 N14-N26
      2 1825 1415 (N13W73 L=234)   739    152        0 (1797/1396-1854/1432)  -4.4 W70-W77 N11-N15
      3 1332 1326 (N19W25 L=186)   413    158        2 (1271/1279-1391/1346)   7.1 W22-W28 N15-N22
      4  347 1269 (N30E48 L=113)   299     63        0 ( 316/1255- 374/1283)   9.5 E52-E45 N29-N31
      5  735 1134 (N17E17 L=145)   467    112        2 ( 708/1072- 765/1190) -90.0 E18-E15 N13-N21
      6 1231  904 (S04W11 L=172)   521    156        1 (1155/ 880-1297/ 935) -12.5 W06-W14 S07-S01
      7 1325  765 (S14W15 L=176)   836    148        4 (1289/ 680-1359/ 879)  71.6 W14-W17 S20-S06
      8 1670  767 (S21W42 L=203)   449    113        1 (1640/ 739-1701/ 795) -26.6 W39-W45 S24-S18
      9  898  554 (S21E16 L=145)   487    204        0 ( 882/ 457- 916/ 639)  63.4 E15-E18 S28-S15
     10 1180  502 (S29W02 L=163)   887    237        2 (1095/ 438-1246/ 583)  30.3 E03-W06 S35-S22
     11  448  450 (S23E56 L=105)   761    186        2 ( 399/ 422- 489/ 486)  12.1 E59-E54 S25-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