Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  808 1787 (N61W17 L= 18)   615     61        2 ( 777/1756- 840/1818) -26.6 W13-W21 N58-N64
      2 1216 1771 (N43W45 L= 47)   786    155        0 (1135/1741-1284/1791) -48.8 W40-W49 N39-N47
      3 1548 1740 (N29W70 L= 72)   726    110        0 (1495/1728-1603/1751) -14.9 W63-W78 N26-N30
      4  628 1667 (N58E09 L=352)   870    133        2 ( 559/1626- 676/1714)  10.6 E18-E02 N55-N61
      5 1314 1598 (N29W38 L= 39)   398     54        0 (1282/1577-1344/1620) -36.9 W35-W39 N27-N31
      6  414 1536 (N55E37 L=324)   577     75        1 ( 391/1499- 436/1571)  51.3 E41-E34 N52-N57
      7 1566 1542 (N18W52 L= 54)   494    143        0 (1525/1501-1602/1566) -71.6 W48-W55 N14-N20
      8 1437 1415 (N16W36 L= 38)   568    254        0 (1314/1371-1541/1446) -45.0 W29-W41 N10-N21
      9  417 1167 (N29E36 L=325)   264     82        0 ( 375/1152- 455/1186)   0.0 E39-E32 N28-N30
     10 1600  687 (S30W27 L= 29)   346     68        0 (1562/ 674-1636/ 705) -18.4 W24-W30 S31-S28
     11 1706  648 (S37W37 L= 38)   319     50        0 (1684/ 629-1729/ 670)   0.0 W34-W39 S38-S35
     12 1053  574 (S20E11 L=351)   277    108        0 (1041/ 519-1070/ 633)  74.1 E12-E08 S23-S16
     13  905  501 (S21E22 L=339)   259    139        0 ( 890/ 437- 921/ 561)  49.4 E25-E19 S24-S17
     14 1417  416 (S43W08 L=  9)  1374    278        2 (1309/ 368-1535/ 491)  -2.9 E01-W19 S45-S40
     15 1027  397 (S31E19 L=342)   271     93        0 ( 997/ 352-1055/ 428)  29.7 E23-E16 S33-S29
     16  630  288 (S28E54 L=308)   912    103        1 ( 577/ 257- 683/ 323) -56.3 E57-E50 S31-S25
     17 1050  244 (S44E27 L=334)   711    104        0 (1025/ 202-1091/ 292)  68.2 E29-E24 S48-S38
  

    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