Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1333 1588 (N41W11 L=192)   494     86        1 (1302/1556-1364/1625) -31.0 W08-W14 N39-N43
      2  324 1532 (N18E67 L=114)   608     43        1 ( 306/1506- 344/1563) -45.0 E70-E63 N15-N19
      3  390 1382 (N10E50 L=131)   319     84        1 ( 368/1356- 425/1407)  45.0 E52-E48 N08-N11
      4 1772 1262 (N28W52 L=233)   421     85        0 (1747/1238-1805/1292)  -9.5 W49-W55 N26-N29
      5 1732 1217 (N24W47 L=228)   702    129        1 (1703/1166-1765/1273) -50.2 W44-W50 N21-N27
      6 1829 1163 (N22W58 L=239)   852    114        2 (1805/1113-1849/1215)  81.9 W55-W61 N19-N26
      7  192 1035 (S14E63 L=118)  3059    374        7 ( 123/ 868- 254/1163) -46.7 E77-E54 S24-S07
      8  712  762 (S22E15 L=166)  1757    458        3 ( 542/ 644- 872/ 889)  47.7 E26-E06 S33-S11
      9 1583  581 (S17W48 L=229)   296     21        0 (1562/ 563-1606/ 598)  18.4 W46-W49 S18-S16
     10 1557  506 (S22W50 L=231)   655    173        2 (1486/ 459-1618/ 552)  48.8 W44-W53 S26-S18
     11 1101  201 (S56W36 L=217)  1547    226        3 ( 992/ 181-1189/ 228)  21.3 W25-W46 S60-S51
  

    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