Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1356 1798 (N49W58 L=328)   829    102        1 (1326/1759-1402/1830)   0.0 W49-W66 N47-N51
      2 1487 1734 (N41W63 L=333)  1558    149        3 (1433/1686-1534/1783)  66.0 W55-W67 N36-N45
      3 1476 1626 (N34W50 L=320)   687    187        0 (1416/1549-1508/1680)  73.3 W46-W52 N29-N39
      4 1425 1599 (N34W44 L=314)   383    120        0 (1389/1555-1464/1643) -90.0 W42-W44 N30-N37
      5 1688 1550 (N24W66 L=336)   581     89        1 (1650/1525-1733/1576) -35.5 W62-W71 N21-N26
      6  579 1486 (N43E27 L=243)   371     89        0 ( 541/1469- 627/1502) -14.0 E30-E22 N41-N44
      7 1736 1350 (N12W58 L=328)   473     69        0 (1718/1323-1758/1393)  45.0 W55-W59 N10-N14
      8 1074 1217 (N17W06 L=277)   904    372        3 ( 910/1156-1234/1272) -31.4 E03-W15 N11-N22
      9  553 1147 (N21E29 L=241)   218     80        1 ( 540/1108- 572/1182)  63.4 E29-E27 N18-N22
     10 1601 1042 (S03W37 L=307)   692    212        2 (1560/ 944-1640/1130)  41.6 W32-W41 S07-N02
     11  907 1081 (N12E06 L=264)   236    116        0 ( 850/1057- 953/1098) -33.7 E09-E03 N09-N13
     12 1866  938 (S16W62 L=332)  1754    302        4 (1768/ 804-1926/1033)  31.4 W49-W73 S22-S11
     13 1001  958 (N03E02 L=268)   227    125        0 ( 980/ 904-1015/1011)  63.4 E04-E01 N00-N06
     14 1902  866 (S22W69 L=339)   781     63        1 (1882/ 839-1924/ 889)  15.9 W64-W74 S23-S20
     15 1567  479 (S40W33 L=303)   328     39        1 (1550/ 451-1588/ 500) -76.0 W30-W34 S42-S38
     16  719  445 (S26E32 L=238)   280     75        1 ( 688/ 417- 746/ 473)  23.2 E35-E28 S27-S23
     17 1457  362 (S48W23 L=294)   538     78        1 (1429/ 330-1492/ 393)  38.7 W20-W26 S50-S45
  

    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