Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1143 1842 (N65W23 L=343)  1853    233        3 (1019/1813-1254/1868)   2.9 W04-W45 N61-N67
      2  427 1611 (N44E58 L=262)   431     79        1 ( 407/1573- 445/1638) -22.6 E64-E52 N40-N45
      3 1431 1548 (N36W35 L=355)   334     75        0 (1391/1531-1467/1563) -18.4 W31-W37 N35-N37
      4  747 1364 (N26E18 L=302)  1170    390        4 ( 683/1269- 838/1458) -61.7 E23-E11 N19-N32
      5 1373 1380 (N25W26 L=346)   413     98        1 (1346/1328-1388/1421)  71.6 W23-W26 N21-N28
      6  369 1352 (N25E50 L=271)   589    111        1 ( 347/1306- 387/1410) -74.1 E52-E47 N22-N29
      7 1908 1155 (N07W75 L= 35)   808     72        0 (1893/1115-1920/1195)  35.5 W70-W79 N04-N09
      8  364 1159 (N13E46 L=274)   423    117        2 ( 316/1130- 412/1189) -11.3 E51-E41 N11-N14
      9  943 1119 (N10E05 L=315)   193     59        2 ( 918/1103- 978/1141) -14.0 E06-E02 N09-N11
     10 1748  859 (S10W51 L= 11)   475     32        0 (1726/ 830-1772/ 881)  56.3 W49-W53 S11-S08
     11  557  798 (S09E31 L=289)  1200    260        2 ( 460/ 742- 636/ 835)  -4.4 E39-E26 S12-S06
     12  409  792 (S10E43 L=277)   457     83        2 ( 389/ 744- 433/ 829)  68.2 E44-E40 S12-S07
     13 1073  679 (S18W02 L=322)   396    146        0 (1056/ 612-1089/ 753)  78.7 W01-W03 S22-S12
     14 1559  646 (S22W37 L=357)   326    118        0 (1536/ 602-1602/ 701)  66.8 W35-W39 S25-S18
     15 1436  540 (S29W29 L=349)   328     63        1 (1417/ 516-1458/ 574)  90.0 W27-W30 S30-S26
     16 1396  477 (S34W27 L=347)   726    135        2 (1354/ 444-1446/ 515) -24.0 W22-W31 S36-S30
  

    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