Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1476 1671 (N28W55 L=239)   971    211        0 (1377/1636-1568/1699) -45.0 W48-W60 N22-N32
      2 1132 1667 (N39W31 L=214)   329     58        0 (1102/1649-1168/1685)   0.0 W27-W34 N38-N40
      3 1654 1537 (N15W59 L=243)   399     62        0 (1619/1519-1687/1554) -26.6 W56-W62 N12-N16
      4 1230 1435 (N23W24 L=208)   432    101        0 (1173/1411-1270/1454) -15.9 W20-W27 N20-N24
      5  642 1074 (N18E21 L=163)   532    192        0 ( 563/1043- 749/1100) -13.0 E26-E13 N16-N19
      6 1862  985 (S22W57 L=241)   412     62        1 (1833/ 963-1896/1001) -18.4 W53-W62 S22-S19
      7 1209  899 (S06W07 L=190)   452    243        1 (1123/ 838-1283/ 942) -56.3 W02-W10 S10-S01
      8 1621  843 (S21W31 L=215)   495    120        0 (1597/ 782-1646/ 899)  45.0 W28-W34 S24-S18
      9  173  781 (N11E70 L=114)  2205    166        4 ( 132/ 715- 204/ 836)  48.8 E79-E65 N06-N14
     10  461  807 (N07E40 L=144)   391     64        2 ( 425/ 789- 500/ 825) -21.8 E42-E37 N06-N08
     11  602  335 (S24E52 L=132)  1256    178        2 ( 564/ 253- 641/ 430)  52.4 E57-E46 S30-S17
     12 1485  359 (S52W13 L=196)  1299    168        1 (1411/ 333-1568/ 396)  -6.7 W05-W22 S54-S48
  

    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