Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1175 1740 (N45W09 L= 98)   297     25        1 (1157/1719-1195/1761)  76.0 W07-W10 N42-N46
      2 1346 1718 (N44W24 L=113)   348     44        1 (1324/1696-1370/1740)  11.3 W21-W26 N42-N46
      3  618 1611 (N32E33 L= 56)   269     83        0 ( 592/1570- 644/1650) -40.6 E36-E29 N29-N35
      4 1379 1508 (N28W23 L=111)   469    160        1 (1333/1458-1430/1578) -53.1 W19-W25 N24-N32
      5 1838 1172 (N10W58 L=147)  4484    726       10 (1779/ 939-1902/1453) -84.3 W54-W66 S04-N27
      6 1742 1422 (N25W53 L=142)   416     56        0 (1716/1392-1769/1451) -31.0 W50-W55 N23-N26
      7  253 1284 (N11E57 L= 32)   558     89        0 ( 239/1241- 269/1330) -71.6 E58-E54 N07-N13
      8 1778 1222 (N12W52 L=141)   533     91        3 (1741/1186-1807/1256) -18.4 W48-W54 N10-N14
      9 1725 1139 (N07W47 L=135)   561    111        2 (1701/1084-1747/1192)  66.8 W44-W48 N03-N10
     10 1806 1097 (N05W54 L=143)   455     96        1 (1780/1070-1844/1128)  23.2 W51-W58 N03-N07
     11  800  877 (S14E13 L= 76)   399    167        1 ( 753/ 815- 863/ 915)  40.6 E16-E09 S17-S11
     12  142  837 (S16E71 L= 18)   789     94        1 ( 103/ 804- 184/ 868)   9.0 E82-E63 S17-S14
     13 1006  776 (S19E00 L= 89)   416     74        2 ( 983/ 740-1035/ 821)  78.7 E01-W01 S21-S16
     14  961  733 (S22E03 L= 86)   195     59        1 ( 926/ 713- 991/ 750) -11.3 E05-W00 S23-S20
     15  840  667 (S27E11 L= 78)   248    101        0 ( 813/ 632- 875/ 713) -51.3 E12-E08 S28-S23
     16 1152  534 (S34W11 L=100)   352     81        1 (1118/ 515-1188/ 558) -21.8 W08-W13 S35-S32
  

    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