Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1371 1776 (N53W49 L= 52)   891     72        0 (1329/1747-1409/1805)  10.0 W40-W57 N51-N55
      2 1260 1767 (N54W34 L= 38)   536     53        2 (1231/1737-1287/1791)  24.0 W30-W39 N51-N56
      3 1602 1670 (N41W66 L= 69)   635     72        1 (1572/1643-1644/1695) -59.0 W61-W69 N38-N43
      4  369 1536 (N41E59 L=305)   563     87        0 ( 324/1522- 411/1553)   4.8 E66-E53 N39-N42
      5  854 1454 (N33E09 L=354)  1169    251        2 ( 760/1397- 946/1539) -35.5 E17-E02 N29-N39
      6  678 1480 (N36E23 L=340)   367    119        0 ( 656/1425- 705/1533)  66.8 E24-E21 N32-N39
      7  883 1217 (N18E08 L=356)   452    200        0 ( 832/1150- 946/1272)  41.2 E11-E03 N13-N20
      8 1354 1052 (N04W21 L= 24)   254     59        1 (1331/1027-1389/1069)  26.6 W19-W23 N02-N05
      9 1738  850 (S12W49 L= 53)   685    109        0 (1703/ 807-1781/ 884)  33.7 W46-W53 S14-S10
     10  382  534 (S24E54 L=309)  1051    199        1 ( 326/ 488- 438/ 570) -35.5 E62-E48 S27-S21
     11  907  507 (S28E12 L=352)   418    104        0 ( 859/ 493- 967/ 522)   7.1 E15-E07 S28-S26
     12 1225  383 (S40W11 L= 14)   286     56        1 (1195/ 366-1258/ 406)  31.0 W08-W13 S41-S38
  

    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