Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  249 1459 (N04E71 L=266)   997    171        3 ( 215/1436- 273/1506)  45.0 E76-E67 N02-N07
      2 1710 1122 (N20W41 L= 18)   298    111        0 (1693/1063-1727/1175) -56.3 W39-W43 N16-N22
      3 1860 1143 (N26W56 L= 32)   510     69        1 (1837/1113-1883/1169) -56.3 W52-W58 N24-N27
      4 1913 1000 (N19W65 L= 42)   659    176        0 (1876/ 989-1946/1014)  16.7 W60-W70 N17-N20
      5  192  910 (S29E57 L=279)   551     67        1 ( 175/ 880- 217/ 935) -56.3 E59-E54 S31-S27
      6  510  891 (S24E29 L=308)   342     93        2 ( 466/ 857- 545/ 921) -15.9 E32-E25 S24-S22
      7 1812  813 (N06W59 L= 36)   557    115        0 (1780/ 759-1841/ 863)  76.0 W57-W61 N01-N09
      8 1592  563 (S14W50 L= 27)   688    230        1 (1491/ 522-1677/ 615)  39.3 W43-W55 S18-S09
      9 1775  529 (S11W72 L= 48)   899    111        0 (1738/ 477-1804/ 579)  74.1 W70-W74 S14-S07
     10 1700  504 (S14W64 L= 40)  1211    235        4 (1644/ 468-1737/ 542)  45.0 W57-W68 S17-S11
     11 1125  385 (S39W28 L=  4)   841    257        2 (1098/ 313-1161/ 482) -37.6 W21-W34 S43-S32
     12 1287  313 (S38W43 L= 20)   998    132        2 (1227/ 289-1339/ 347)  20.6 W38-W47 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