Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  921 1903 (N48E51 L=220)   731     83        1 ( 899/1872- 950/1922) -40.6 E56-E45 N44-N50
      2  792 1824 (N36E48 L=223)   618     83        2 ( 772/1781- 812/1864) -41.6 E53-E43 N32-N40
      3 1097 1628 (N31E15 L=256)   536    108        1 (1068/1571-1116/1681) -63.4 E17-E12 N27-N35
      4 1349 1551 (N35W04 L=275)  1014    302        1 (1271/1423-1404/1653) -35.5 E04-W10 N27-N40
      5 1483 1175 (N16W22 L=293)   824    201        3 (1428/1098-1526/1278) -41.6 W17-W26 N12-N21
      6  492 1156 (S12E35 L=236)   447     94        1 ( 473/1104- 514/1208) -63.4 E37-E33 S15-S09
      7  748 1004 (S15E15 L=256)   234    140        0 ( 716/ 953- 796/1044)  74.1 E15-E12 S18-S11
      8 1452  977 (N04W25 L=296)   345    108        0 (1409/ 936-1480/1024)  74.1 W23-W26 N00-N07
      9  647  960 (S20E20 L=251)   216     61        0 ( 627/ 929- 674/ 989) -26.6 E21-E17 S20-S18
     10  988  410 (S42W21 L=292)   578    200        0 ( 889/ 392-1084/ 433)  18.4 W12-W27 S45-S39
     11 1309  403 (S31W42 L=313)   587    113        0 (1254/ 387-1363/ 429)  12.5 W36-W45 S33-S29
  

    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