Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  804 1644 (N40E26 L=124)   765    113        2 ( 753/1602- 853/1679) -16.7 E31-E21 N36-N42
      2 1429 1593 (N43W28 L=178)   949    157        3 (1374/1543-1486/1643) -33.7 W23-W32 N39-N47
      3  374 1515 (N25E58 L= 91)  1638    191        4 ( 347/1428- 401/1601) -32.9 E70-E52 N19-N30
      4 1649 1499 (N38W48 L=198)   557    110        0 (1612/1461-1681/1529) -51.3 W44-W51 N35-N40
      5  851 1389 (N22E15 L=135)   544    225        1 ( 821/1287- 890/1487)  85.9 E16-E14 N15-N29
      6 1679 1225 (N20W45 L=195)   283    136        0 (1650/1163-1708/1279) -60.3 W42-W47 N16-N23
      7  500  858 (S14E33 L=117)  1989    443        3 ( 380/ 686- 621/1000)  49.2 E44-E24 S26-S04
      8 1492  711 (S14W34 L=184)   454    136        0 (1428/ 688-1553/ 731)   0.0 W29-W39 S14-S12
      9 1750  670 (S14W59 L=209)   492     64        2 (1722/ 637-1783/ 693)  -7.1 W54-W62 S16-S12
     10  329  644 (S30E52 L= 98)   392     55        1 ( 298/ 630- 356/ 665)   7.1 E56-E48 S31-S28
     11 1413  507 (S28W34 L=184)   649    202        0 (1324/ 462-1485/ 558)  45.0 W28-W38 S31-S23
     12  584  398 (S46E32 L=118)   556     76        1 ( 547/ 377- 625/ 424)  24.4 E37-E26 S48-S43
  

    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