Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1602 1458 (N31W45 L=334)  1236    216        0 (1525/1374-1676/1526) -50.7 W39-W50 N24-N35
      2  382 1488 (N31E54 L=235)   365     58        1 ( 359/1468- 414/1514) -14.0 E57-E49 N29-N32
      3 1755 1469 (N31W64 L=353)   545     80        0 (1720/1440-1792/1491) -20.6 W60-W68 N28-N32
      4 1680  961 (S01W45 L=334)  1213    184        5 (1626/ 895-1751/1010)  32.0 W40-W51 S04-N01
      5 1816  840 (S09W59 L=348)   627     85        1 (1794/ 783-1834/ 880) -71.6 W56-W62 S12-S06
      6  254  770 (S15E58 L=231)   534    136        2 ( 229/ 730- 310/ 815) -20.6 E60-E52 S17-S11
      7 1161  689 (S18W09 L=298)   235     64        0 (1141/ 656-1178/ 724) -63.4 W07-W10 S19-S15
      8  191  696 (S20E71 L=218)   978     66        0 ( 163/ 670- 218/ 721)  13.2 E82-E64 S22-S18
      9  315  630 (S24E55 L=234)   521    122        0 ( 288/ 576- 340/ 680)  37.9 E60-E50 S27-S20
     10 1427  555 (S27W29 L=318)   323    133        0 (1404/ 489-1442/ 597) -63.4 W26-W30 S31-S23
     11 1749  547 (S29W63 L=352)   779     94        0 (1704/ 519-1786/ 587)  26.6 W57-W67 S30-S25
     12 1092  516 (S29W05 L=294)  1429    288        3 ( 962/ 472-1218/ 577)  23.6 E03-W13 S32-S25
  

    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