Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  684 1710 (N57E22 L= 75)   648     99        1 ( 628/1689- 732/1731)  -4.4 E30-E16 N55-N59
      2  475 1695 (N58E55 L= 43)  1003    166        2 ( 406/1657- 545/1715)   7.8 E67-E43 N55-N59
      3  520 1364 (N33E32 L= 65)  1089    325        1 ( 435/1235- 586/1465) -50.5 E42-E26 N23-N40
      4 1803 1368 (N14W65 L=163)   485     56        0 (1778/1348-1830/1385) -31.0 W62-W68 N12-N15
      5  389 1282 (N28E43 L= 54)   287     87        0 ( 372/1243- 416/1325) -50.2 E45-E40 N25-N31
      6 1750 1229 (N07W53 L=151)   506    107        0 (1723/1185-1772/1279) -80.5 W51-W55 N04-N10
      7  989 1111 (N12E01 L= 96)   202     54        2 ( 965/1082-1010/1129) -56.3 E02-W00 N09-N13
      8 1569  819 (S14W32 L=130)   481    108        0 (1541/ 768-1593/ 868)  51.3 W30-W34 S17-S11
      9 1465  555 (S30W24 L=122)  2244    436        2 (1280/ 472-1646/ 650)  11.7 W10-W39 S33-S26
     10  567  604 (S15E36 L= 62)   646    136        0 ( 506/ 572- 635/ 639)  10.3 E41-E30 S16-S12
     11 1130  365 (S39E02 L= 95)   442    152        0 (1053/ 339-1192/ 394)  11.3 E08-W02 S39-S37
     12  701  206 (S48E49 L= 48)  1056     77        2 ( 663/ 185- 747/ 233)   7.6 E58-E42 S49-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