Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  467 1669 (N34E49 L=132)  1140    259        0 ( 420/1627- 522/1723) -28.1 E57-E42 N30-N38
      2 1609 1417 (N24W37 L=218)   701    159        1 (1536/1390-1668/1438)   0.0 W31-W42 N22-N26
      3 1821 1289 (N19W55 L=236)   530     63        1 (1796/1258-1845/1327) -53.1 W52-W57 N17-N21
      4  387 1265 (N06E44 L=137)   659     85        2 ( 366/1218- 417/1313) -63.4 E45-E40 N03-N09
      5 1916 1124 (N10W66 L=248)  1570    313        4 (1874/1030-1952/1189) -41.6 W60-W73 N04-N14
      6 1294 1000 (S04W16 L=197)  1070    322        1 (1238/ 859-1331/1149) -76.0 W11-W18 S12-N04
      7 1669  952 (S04W42 L=224)   401    117        0 (1616/ 933-1712/ 970)  15.9 W38-W46 S05-S02
      8  789  883 (S15E13 L=168)   195     94        0 ( 745/ 859- 844/ 897)  18.4 E16-E10 S16-S13
      9 1661  802 (S12W44 L=225)   385     99        1 (1624/ 780-1700/ 818)  18.4 W41-W47 S14-S11
     10  439  791 (S22E38 L=143)   293     93        0 ( 405/ 758- 482/ 816)  26.6 E40-E34 S23-S20
     11 1811  645 (S19W63 L=244)  2959    365        4 (1702/ 568-1899/ 714)  34.5 W54-W73 S25-S14
     12 1460  418 (S39W40 L=221)   806    207        2 (1375/ 397-1554/ 449)  21.0 W33-W46 S40-S35
  

    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