Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1284 1540 (N34W03 L= 77)  1482    250        3 (1223/1434-1334/1641) -45.0 E04-W08 N28-N40
      2  383 1527 (N12E59 L= 15)   379     50        0 ( 367/1495- 399/1552) -53.1 E61-E56 N09-N13
      3  668 1459 (N14E33 L= 41)   285     98        0 ( 649/1416- 700/1495) -80.5 E34-E31 N10-N16
      4  368 1467 (N08E57 L= 17)   460     65        1 ( 347/1443- 402/1494)   0.0 E59-E53 N06-N09
      5 1591 1168 (N19W32 L=106)  1115    215        2 (1541/1067-1635/1260)  85.6 W29-W34 N12-N25
      6 1819 1036 (N18W55 L=129)  1515    178        2 (1782/ 946-1848/1120) -56.3 W50-W59 N12-N22
      7 1292 1044 (N04W15 L= 89)   327     78        0 (1251/1023-1325/1070)  45.0 W12-W16 N02-N06
      8 1900  868 (N11W71 L=144)   694     71        1 (1870/ 844-1921/ 901)  51.3 W65-W75 N08-N13
      9  938  743 (S21W02 L= 76)   452    238        0 ( 843/ 706-1053/ 788)  42.0 E03-W07 S25-S15
     10  542  714 (S33E24 L= 49)  1030    250        1 ( 446/ 679- 634/ 768)  23.2 E32-E18 S35-S27
     11 1556  559 (S16W48 L=122)   841    212        3 (1479/ 513-1626/ 604)  48.8 W43-W52 S20-S12
     12 1655  553 (S14W57 L=131)   358     54        1 (1625/ 531-1680/ 577)  45.0 W54-W59 S15-S11
  

    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