Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  778 1817 (N57E40 L=308)  1376    114        2 ( 740/1751- 821/1853) -54.0 E47-E33 N50-N61
      2  665 1822 (N55E57 L=291)   618     37        0 ( 635/1805- 692/1836)  -4.1 E64-E50 N54-N56
      3  960 1596 (N39E10 L=338)   293     60        0 ( 939/1571- 988/1613) -14.0 E11-E07 N37-N40
      4  514 1494 (N28E43 L=305)   684    136        1 ( 448/1462- 570/1537) -13.0 E50-E37 N26-N30
      5 1503 1408 (N29W32 L= 20)   671    211        0 (1447/1337-1572/1480) -41.6 W28-W37 N25-N34
      6  319 1397 (N20E58 L=290)  1526    178        3 ( 264/1324- 361/1475) -24.0 E68-E50 N15-N24
      7  778 1323 (N19E18 L=330)   548    197        0 ( 696/1260- 844/1363) -24.4 E24-E13 N15-N21
      8 1345 1077 (N07W20 L=  8)   436     77        1 (1315/1038-1373/1110)  53.1 W18-W21 N04-N09
      9 1260  699 (S17W17 L=  5)   293     82        0 (1232/ 662-1284/ 734)  68.2 W16-W18 S19-S14
     10 1449  358 (S40W44 L= 32)   359     46        0 (1433/ 334-1467/ 387) -59.0 W40-W46 S42-S37
     11 1098  318 (S46W13 L=  1)   343     58        2 (1070/ 298-1125/ 338) -31.0 W10-W15 S47-S44
     12 1258  298 (S47W29 L= 17)   456     16        0 (1236/ 277-1284/ 320) -11.3 W26-W31 S48-S44
  

    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