Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1465 1553 (N32W42 L=319)  1535    349        1 (1309/1494-1608/1592) -25.5 W30-W52 N25-N36
      2  483 1534 (N45E40 L=237)   327     72        0 ( 455/1501- 507/1564) -32.0 E45-E36 N42-N47
      3  330 1367 (N34E52 L=225)  1825    305        1 ( 216/1285- 428/1429) -15.3 E73-E40 N28-N38
      4  937 1342 (N26E02 L=275)   484    204        0 ( 841/1330-1035/1352) -14.0 E08-W04 N24-N27
      5 1577 1235 (N11W39 L=316)   798    239        3 (1543/1123-1624/1330) -90.0 W37-W41 N03-N16
      6  797 1163 (N17E13 L=264)   193     67        0 ( 771/1140- 830/1180)  26.6 E14-E10 N15-N17
      7 1241 1145 (N10W15 L=292)   248     75        0 (1208/1119-1273/1162) -45.0 W12-W16 N08-N11
      8 1828 1106 (S02W59 L=336)  1243    196        4 (1777/1055-1893/1155) -12.1 W54-W69 S04-N02
      9  802 1037 (N09E13 L=264)   340    146        0 ( 788/ 976- 822/1111) -90.0 E14-E12 N05-N13
     10 1731  845 (S16W47 L=324)  1999    364        5 (1692/ 710-1779/ 964)  86.6 W43-W52 S24-S07
     11 1831  935 (S12W58 L=334)   597    110        1 (1776/ 902-1866/ 962) -24.4 W51-W62 S14-S09
     12 1662  872 (S13W41 L=317)   680    164        2 (1628/ 807-1693/ 950)  66.0 W37-W43 S16-S07
     13  361  736 (S06E49 L=228)  1528    314        3 ( 306/ 654- 413/ 822) -78.7 E54-E45 S10-N00
     14 1282  744 (S15W13 L=289)   329     94        0 (1238/ 729-1331/ 761) -18.4 W09-W15 S16-S13
     15  299  619 (S13E60 L=217)  4774    406       15 ( 186/ 451- 389/ 754) -84.6 E72-E53 S25-S04
     16 1192  614 (S22W06 L=282)  1313    317        2 (1155/ 507-1265/ 730)  86.2 W02-W09 S29-S14
     17  958  441 (S31E13 L=264)   338     75        1 ( 918/ 415- 985/ 474)  33.7 E16-E10 S32-S28
     18 1381  378 (S43W19 L=296)   840    112        3 (1322/ 347-1427/ 413)   0.0 W13-W23 S45-S40
  

    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