Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1243 1695 (N41W10 L=342)   304     93        0 (1204/1669-1291/1713) -15.9 W07-W14 N39-N42
      2 1506 1603 (N37W32 L=  3)  1608    250        5 (1394/1532-1585/1653) -11.3 W22-W37 N32-N40
      3  415 1570 (N26E52 L=280)   540    100        2 ( 361/1548- 454/1586)  -5.2 E58-E47 N24-N26
      4  582 1545 (N24E36 L=296)   400     77        1 ( 552/1508- 607/1569) -68.2 E38-E33 N21-N26
      5  730 1514 (N23E24 L=308)   311    120        0 ( 687/1467- 778/1542)  45.0 E26-E20 N19-N25
      6  425 1520 (N22E48 L=283)   492     64        0 ( 388/1505- 469/1540)   0.0 E52-E43 N21-N23
      7 1198 1371 (N18W08 L=339)   503    212        0 (1096/1345-1291/1408)  -4.8 W01-W13 N16-N19
      8  968 1358 (N15E06 L=325)   572    200        0 ( 898/1299-1025/1388) -26.6 E10-E02 N11-N16
      9 1835 1298 (N21W59 L= 31)  2245    203        4 (1778/1213-1875/1381) -63.4 W53-W64 N15-N25
     10  734 1305 (N10E20 L=311)   622    128        0 ( 694/1257- 776/1364) -45.0 E23-E16 N07-N13
     11 1240 1245 (N10W11 L=343)   213     48        0 (1219/1224-1266/1271) -33.7 W09-W12 N08-N11
     12 1080 1064 (S02W03 L=334)   839    216        2 (1034/ 967-1145/1146)  70.0 W00-W06 S08-N03
     13 1408 1002 (S03W24 L=355)   277     71        1 (1387/ 969-1444/1027)  45.0 W22-W25 S04-S01
     14 1613  734 (S16W43 L= 14)   261     42        0 (1586/ 721-1648/ 747)   0.0 W40-W45 S17-S15
     15  721  578 (S35E18 L=314)  1944    401        3 ( 566/ 490- 859/ 721)  37.7 E32-E07 S42-S24
     16  370  609 (S35E49 L=283)   403     72        2 ( 339/ 580- 408/ 632)  18.4 E53-E44 S36-S33
     17 1247  291 (S51W32 L=  3)   362     61        1 (1211/ 275-1280/ 308)  31.0 W28-W34 S52-S49
  

    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