Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1265 1588 (N39W05 L=149)   255     55        0 (1243/1568-1297/1603)  26.6 W03-W07 N37-N40
      2  928 1474 (N24E15 L=129)   624    243        0 ( 849/1384-1015/1553) -28.1 E22-E07 N19-N27
      3 1379 1530 (N38W14 L=158)   303     60        0 (1359/1507-1411/1549)  45.0 W12-W16 N35-N39
      4 1674 1380 (N34W40 L=184)  2125    388        3 (1540/1310-1798/1466)  -7.4 W28-W52 N30-N37
      5 1399 1217 (N17W19 L=164)   773    233        1 (1345/1141-1488/1258)  36.9 W15-W25 N12-N20
      6 1383 1084 (N09W20 L=165)   298    114        1 (1350/1049-1438/1132) -15.9 W17-W24 N07-N11
      7 1135  983 (S02W07 L=151)   445    155        0 (1041/ 955-1186/1026)  35.5 W02-W09 S04-N01
      8  493  710 (S31E30 L=114)   434    139        0 ( 460/ 654- 518/ 762)  66.0 E32-E28 S35-S26
      9 1145  697 (S18W14 L=158)   246     39        0 (1124/ 676-1166/ 718) -33.7 W12-W15 S19-S17
     10 1606  664 (S11W46 L=190)   823    174        1 (1545/ 637-1671/ 700)   5.2 W40-W51 S12-S09
     11 1495  382 (S31W51 L=195)   709    155        4 (1457/ 362-1545/ 409)   6.3 W46-W55 S32-S29
     12 1238  253 (S47W41 L=185)   372     50        0 (1214/ 234-1256/ 277)  90.0 W38-W42 S48-S44
  

    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