Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  677 1773 (N33E47 L=360)   656     88        0 ( 642/1743- 717/1799) -16.7 E51-E41 N31-N34
      2 1255 1754 (N46E03 L= 44)   511     82        1 (1220/1724-1287/1789)  63.4 E04-W00 N42-N49
      3  455 1668 (N22E59 L=348)  3959    243        7 ( 353/1617- 535/1757)  71.6 E72-E50 N17-N29
      4  869 1658 (N28E26 L= 21)   720    233        1 ( 779/1586- 966/1698) -14.0 E33-E17 N25-N29
      5  749 1588 (N21E30 L= 16)   259    164        0 ( 689/1556- 809/1612)   7.1 E33-E25 N18-N22
      6 1520 1465 (N31W23 L= 69)  2145    493        3 (1331/1304-1657/1592) -19.7 W06-W34 N24-N35
      7 1112 1426 (N18E04 L= 43)   538    212        0 (1038/1355-1180/1505) -31.0 E09-W02 N15-N21
      8  678 1243 (S01E25 L= 22)   352     81        0 ( 649/1212- 700/1290)  76.0 E26-E23 S03-N01
      9 1809 1142 (N19W50 L= 97)  2301    292       13 (1739/1052-1901/1263) -19.3 W42-W62 N14-N24
     10 1838  966 (N10W56 L=103)  1224    162        2 (1813/ 882-1860/1048) -84.3 W53-W58 N04-N14
     11  509  931 (S22E31 L= 16)   301     80        1 ( 498/ 884- 523/ 973) -68.2 E31-E29 S24-S19
     12  162  827 (S33E68 L=338)  1246    185        3 ( 133/ 745- 211/ 901)  26.6 E79-E58 S38-S28
     13 1714  650 (S11W56 L=102)   395     56        2 (1699/ 623-1733/ 674)  71.6 W53-W58 S12-S09
     14  530  626 (S41E27 L= 20)   283     76        0 ( 494/ 606- 571/ 640)  23.2 E30-E23 S43-S39
     15 1376  529 (S27W36 L= 82)   623    109        1 (1336/ 491-1429/ 571) -11.3 W31-W41 S28-S25
     16 1486  323 (S35W60 L=106)   936     97        1 (1439/ 291-1530/ 361)   0.0 W51-W65 S37-S33
     17 1386  271 (S42W57 L=104)   590     65        0 (1364/ 233-1406/ 305) -14.0 W51-W63 S43-S39
  

    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