Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1189 1380 (N22W02 L=180)  1059    222        3 (1117/1278-1245/1474) -39.3 E04-W07 N17-N26
      2  202 1002 (S20E58 L=120)   565    161        2 ( 167/ 965- 243/1046)  40.6 E62-E54 S23-S17
      3 1777  891 (N07W54 L=232)  1036    150        1 (1756/ 816-1794/ 965) -69.4 W51-W57 N03-N11
      4  528  862 (S22E28 L=150)   430    165        1 ( 463/ 840- 569/ 917)  45.0 E32-E25 S24-S18
      5 1319  824 (S07W22 L=200)   383    102        1 (1281/ 773-1336/ 855) -45.0 W18-W23 S10-S05
      6 1116  574 (S27W17 L=195)   546    178        0 (1032/ 535-1165/ 643)  53.1 W12-W18 S30-S21
      7 1631  583 (S14W52 L=230)   917    243        1 (1556/ 520-1702/ 628)  55.0 W48-W56 S19-S09
      8 1797  596 (S09W70 L=248)   588    136        3 (1774/ 563-1820/ 624) -23.2 W65-W73 S10-S06
      9  598  600 (S38E20 L=159)   268    145        0 ( 535/ 573- 654/ 618)  33.7 E24-E15 S41-S35
     10 1496  582 (S17W41 L=220)   296     49        1 (1468/ 567-1524/ 602)  33.7 W39-W43 S18-S15
  

    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