Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  493 1716 (N28E58 L=289)  3191    269       13 ( 414/1625- 570/1812) -69.4 E66-E52 N20-N36
      2 1269 1685 (N41W02 L=349)   366     77        0 (1241/1662-1303/1704)   0.0 W00-W04 N38-N42
      3 1412 1655 (N42W14 L=  1)   472     95        1 (1370/1616-1445/1684) -23.2 W10-W17 N39-N43
      4  475 1557 (N16E49 L=299)   778    202        1 ( 419/1480- 535/1632)  79.7 E50-E46 N10-N21
      5  712 1356 (N07E25 L=322)   246     92        0 ( 677/1321- 743/1391)  59.0 E26-E23 N04-N09
      6 1342 1292 (N15W14 L=  1)   572    219        0 (1231/1254-1416/1315)   9.5 W06-W18 N12-N16
      7 1728 1268 (N23W42 L= 30)   302     54        1 (1704/1242-1749/1298) -33.7 W40-W44 N21-N24
      8 1188 1251 (N09W05 L=352)   233     37        1 (1169/1224-1207/1279) -76.0 W04-W06 N07-N11
      9  339 1121 (S12E47 L=300)   609    167        2 ( 272/1065- 391/1165)  37.9 E52-E42 S16-S08
     10 1411 1114 (N06W21 L=  8)   245     95        0 (1365/1091-1434/1146)  45.0 W18-W22 N04-N08
     11 1644 1074 (N09W37 L= 25)   335    122        0 (1633/1018-1674/1115)  81.9 W36-W39 N05-N12
     12  501  993 (S17E32 L=315)   219     96        0 ( 491/ 943- 512/1040)  90.0 E33-E31 S20-S14
     13  380  747 (S34E42 L=305)   817    250        2 ( 282/ 700- 466/ 823)  27.8 E53-E33 S38-S28
     14 1771  720 (S07W58 L= 45)   446     50        1 (1750/ 687-1797/ 745) -76.0 W55-W60 S09-S05
     15 1749  573 (S15W64 L= 51)  4873    328       12 (1606/ 488-1835/ 663)  41.2 W54-W78 S23-S09
     16 1475  543 (S25W42 L= 29)   893    225        1 (1430/ 450-1524/ 653) -72.9 W39-W45 S30-S17
     17 1568  370 (S32W62 L= 50)   598     76        0 (1539/ 347-1602/ 392)  76.0 W58-W65 S33-S29
  

    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