Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  555 1608 (N15E47 L=281)   382     86        0 ( 540/1571- 570/1644) -59.0 E49-E45 N12-N17
      2 1295 1539 (N31W01 L=329)   397     94        1 (1250/1507-1344/1561)   0.0 E02-W04 N29-N32
      3 1698 1511 (N46W34 L=  2)   803    128        2 (1673/1460-1716/1555)  71.6 W31-W36 N41-N50
      4 1527 1424 (N32W20 L=348)   352     85        1 (1502/1393-1557/1469) -31.0 W17-W22 N29-N34
      5  506 1417 (N02E42 L=286)   440     98        2 ( 487/1375- 535/1457) -53.1 E44-E39 N00-N04
      6  703 1379 (N05E28 L=300)   310    100        0 ( 656/1342- 742/1402)   0.0 E31-E24 N03-N05
      7  252 1207 (S15E56 L=272)  2213    466        4 ( 139/1169- 376/1231)   4.4 E70-E44 S18-S11
      8 1823  983 (N16W53 L= 21)   620    106        1 (1802/ 918-1837/1029) -51.3 W49-W54 N11-N18
      9  826  863 (S21E07 L=321)   269     87        0 ( 811/ 821- 840/ 916) -53.1 E08-E05 S22-S18
     10 1762  834 (N06W52 L= 20)   738    159        1 (1718/ 798-1809/ 883)  -5.2 W46-W57 N04-N06
     11 1349  791 (S10W25 L=353)   201     77        0 (1332/ 748-1366/ 834) -78.7 W23-W25 S12-S07
     12 1886  785 (N08W69 L= 37)   784    184        1 (1854/ 762-1922/ 809)  32.0 W63-W74 N05-N10
     13 1508  689 (S10W38 L=  7)   316     73        2 (1478/ 668-1537/ 708)  56.3 W36-W40 S12-S08
     14  826  628 (S34E00 L=328)   202     72        0 ( 787/ 617- 863/ 640)  36.9 E02-W02 S36-S32
     15  992  570 (S33W13 L=341)  1259    272        3 ( 847/ 531-1104/ 630)  40.6 W03-W17 S38-S26
     16 1584  586 (S13W48 L= 16)   363     99        0 (1540/ 571-1627/ 606)  38.7 W45-W51 S15-S10
     17 1434  414 (S27W47 L= 15)  1813    408        3 (1311/ 338-1534/ 486)   2.4 W35-W59 S28-S24
  

    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