Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1483 1637 (N26W53 L=341)   964    232        1 (1399/1564-1558/1721) -90.0 W51-W54 N19-N33
      2  885 1553 (N41W08 L=297)   537     82        1 ( 843/1518- 923/1582)   8.1 W04-W11 N39-N42
      3 1022 1537 (N36W16 L=305)   548    117        2 ( 964/1502-1067/1563) -23.2 W13-W20 N33-N38
      4  333 1516 (N56E52 L=237)  1007    110        1 ( 287/1486- 393/1537)  -3.0 E60-E40 N54-N57
      5 1703 1399 (N06W55 L=344)   913    130        3 (1668/1342-1736/1456)  29.1 W49-W59 N04-N09
      6 1381 1347 (N14W30 L=319)   373     66        1 (1350/1325-1416/1374) -68.2 W28-W31 N11-N16
      7  693 1228 (N26E15 L=274)   241     92        2 ( 659/1194- 723/1251) -36.9 E16-E12 N24-N28
      8 1148 1215 (N14W12 L=301)   373     86        1 (1102/1190-1193/1246) -45.0 W09-W14 N11-N16
      9  196 1175 (N35E62 L=227)   802    106        1 ( 155/1135- 237/1209) -23.2 E69-E55 N31-N37
     10  748 1135 (N20E13 L=275)   220     88        0 ( 715/1094- 769/1164) -80.5 E14-E12 N16-N22
     11 1704 1070 (S11W43 L=332)   357     48        1 (1681/1047-1736/1092)  11.3 W41-W46 S11-S09
     12  726  917 (N08E20 L=269)   445    165        2 ( 648/ 894- 785/ 949) -18.4 E24-E15 N06-N10
     13 1653  804 (S24W33 L=322)   319    131        0 (1610/ 745-1694/ 835)  12.5 W28-W37 S26-S23
     14 1558  642 (S31W23 L=312)   339     75        0 (1542/ 599-1577/ 674)  59.0 W21-W25 S34-S29
     15 1364  577 (S29W08 L=296)   943    214        0 (1289/ 512-1423/ 670)  26.6 W01-W13 S31-S24
     16 1704  579 (S42W37 L=326)  1195    105        3 (1672/ 527-1739/ 631)  74.1 W33-W41 S45-S37
     17 1061  228 (S45E28 L=261)  1657    259        3 ( 943/ 175-1161/ 289)   2.4 E41-E16 S47-S42
  

    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