Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1021 1784 (N40E29 L=343)   823    175        1 ( 939/1734-1090/1821)  -3.6 E37-E21 N38-N42
      2  579 1617 (N16E47 L=325)   675    178        1 ( 554/1538- 610/1695) -55.0 E50-E42 N11-N21
      3  685 1620 (N19E40 L=332)   658    184        2 ( 656/1531- 708/1694) -48.0 E45-E35 N13-N23
      4 1503 1451 (N34W17 L= 29)   901    163        4 (1479/1385-1540/1528) -73.3 W14-W20 N28-N38
      5  779 1410 (N09E25 L=347)   998    282        4 ( 699/1324- 888/1493)  71.6 E27-E21 N01-N16
      6  204 1239 (S14E65 L=307)   720    124        2 ( 155/1216- 262/1262)  12.1 E72-E57 S15-S12
      7  979 1172 (N01E07 L=  5)   658    227        3 ( 897/1143-1062/1225)  16.7 E12-E02 S01-N04
      8 1586 1093 (N14W31 L= 43)   231    148        0 (1549/1028-1621/1148) -35.5 W27-W34 N11-N16
      9  571 1076 (S15E28 L=344)   250     76        0 ( 530/1050- 601/1100)  45.0 E30-E26 S17-S13
     10 1484 1019 (N07W26 L= 38)   392    164        1 (1469/ 936-1501/1089) -49.4 W23-W29 N03-N10
     11  210  874 (S36E59 L=313)  1254    222        3 ( 133/ 798- 286/ 948)  23.2 E77-E49 S41-S29
     12  911  724 (S26W02 L= 14)  1038    351        4 ( 766/ 675-1071/ 792)  47.0 E05-W09 S32-S17
     13 1502  639 (S12W41 L= 53)   319    116        1 (1443/ 618-1555/ 656)  18.4 W36-W45 S14-S11
     14  317  571 (S54E51 L=321)   425     76        0 ( 292/ 555- 355/ 597)  13.0 E56-E43 S55-S52
     15 1362  508 (S24W38 L= 50)   393    146        1 (1294/ 497-1436/ 520)  29.1 W34-W43 S26-S21
     16 1218  343 (S38W41 L= 53)  1919    315        3 (1089/ 269-1351/ 426)  78.1 W36-W44 S47-S28
  

    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