Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1438 1791 (N49W65 L=119)   585     57        0 (1425/1757-1454/1822)  13.2 W58-W75 N46-N51
      2 1761 1484 (N22W69 L=123)   815     87        0 (1716/1456-1802/1511) -29.1 W64-W75 N19-N24
      3  674 1425 (N35E21 L= 33)   414    108        0 ( 637/1371- 698/1468) -60.3 E23-E18 N31-N38
      4 1102 1400 (N28W10 L= 64)   697    188        0 (1025/1343-1161/1430) -24.0 W04-W13 N23-N30
      5 1544 1318 (N16W39 L= 93)   255    100        0 (1534/1270-1563/1367)  56.3 W37-W42 N13-N19
      6  873 1331 (N27E06 L= 48)   307    136        0 ( 811/1312- 935/1350) -20.6 E10-E02 N24-N28
      7  291 1219 (N24E54 L=  1)   405     81        0 ( 277/1178- 311/1261) -56.3 E55-E50 N21-N27
      8  644 1167 (N19E23 L= 31)   221     81        0 ( 606/1149- 681/1187)  21.8 E25-E20 N17-N19
      9 1589  962 (S06W36 L= 90)  1647    352        3 (1525/ 805-1622/1104) -80.5 W32-W38 S15-N03
     10 1829  942 (S11W58 L=112)  2413    385        5 (1725/ 824-1917/1083)  29.7 W46-W72 S16-S04
     11  748  889 (N01E19 L= 35)   350     83        1 ( 708/ 862- 783/ 923)  21.8 E21-E16 N00-N02
     12 1183  721 (S15W06 L= 60)   280    114        0 (1130/ 686-1226/ 743)  15.9 W02-W09 S16-S13
     13  367  526 (S20E58 L=356)  1559    202        3 ( 303/ 439- 402/ 607)  80.5 E61-E54 S26-S14
     14  295  517 (S21E69 L=345)  2204    138        5 ( 250/ 459- 333/ 574) -60.3 E77-E61 S25-S16
     15 1321  493 (S33W14 L= 68)   347     38        0 (1300/ 469-1343/ 518)  33.7 W12-W15 S34-S30
     16 1233  179 (S63W06 L= 60)   581     51        0 (1206/ 160-1262/ 192) -26.6 W02-W10 S66-S60
  

    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