Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1090 1767 (N45E05 L=259)  1374    261        0 ( 951/1728-1194/1806)  -5.4 E18-W03 N43-N47
      2 1558 1628 (N41W36 L=300)   480     83        0 (1514/1604-1596/1646)  -7.1 W31-W39 N40-N42
      3 1363 1504 (N28W17 L=282)   349     71        1 (1325/1479-1395/1546) -53.1 W14-W19 N26-N31
      4  829 1465 (N19E17 L=247)   980    260        0 ( 742/1373- 927/1543) -31.0 E24-E09 N14-N23
      5  253 1497 (N18E68 L=197)   968     86        1 ( 233/1467- 288/1540) -24.4 E73-E60 N16-N21
      6 1294 1458 (N24W13 L=277)   357     84        0 (1256/1430-1332/1491) -26.6 W09-W15 N22-N26
      7 1515 1427 (N25W29 L=293)   365     47        0 (1492/1398-1543/1458) -45.0 W26-W30 N23-N27
      8 1583 1343 (N21W33 L=298)   887    246        1 (1545/1231-1622/1454)  75.1 W30-W37 N13-N28
      9 1708 1219 (N15W44 L=308)   724    162        2 (1665/1150-1750/1293) -45.0 W40-W47 N11-N19
     10  112 1051 (S10E73 L=192)   708    101        0 (  79/1007- 149/1088)  14.0 E82-E66 S11-S07
     11 1202  851 (S14W13 L=277)   243     57        0 (1187/ 817-1221/ 888) -53.1 W11-W14 S15-S11
     12  240  740 (S29E59 L=205)  1016    141        1 ( 208/ 674- 279/ 806)  36.0 E64-E52 S32-S24
     13 1573  745 (S15W39 L=304)   255     76        0 (1530/ 732-1616/ 757)   9.5 W36-W42 S15-S14
     14 1710  722 (S14W52 L=316)   384     53        2 (1689/ 691-1730/ 744)  76.0 W49-W53 S16-S12
     15 1773  682 (S15W60 L=324)  1893    138        4 (1726/ 613-1834/ 735) -69.4 W54-W67 S19-S11
     16 1004  630 (S30W03 L=268)  1285    353        3 ( 884/ 519-1132/ 731)  50.9 E04-W10 S38-S21
  

    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