Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1289 1280 (N14W22 L=278)  1369    462        4 (1173/1188-1465/1311)  -8.1 W12-W33 N09-N17
      2 1831 1023 (S17W53 L=309)  1194    213        2 (1783/ 982-1896/1084) -17.1 W47-W60 S20-S13
      3  620  993 (N15E24 L=231)  1849    449        7 ( 443/ 954- 805/1040) -23.5 E36-E13 N09-N19
      4 1709  826 (S24W38 L=294)   721    163        1 (1664/ 746-1739/ 911)  45.0 W33-W41 S28-S20
      5  319  732 (N06E54 L=202)  1491    297        2 ( 282/ 627- 375/ 887)  68.2 E57-E47 N00-N15
      6 1271  785 (S13W08 L=264)   369    159        1 (1236/ 711-1312/ 863)  41.2 W03-W11 S16-S09
      7 1172  795 (S10W02 L=258)   596    168        2 (1125/ 708-1232/ 851)  35.5 E01-W06 S14-S07
      8 1688  641 (S36W34 L=290)   661    172        2 (1636/ 603-1730/ 666) -37.9 W29-W38 S40-S33
      9 1281  549 (S28W02 L=258)   319     65        0 (1247/ 524-1312/ 581)  26.6 W00-W04 S29-S26
     10  698  407 (S21E39 L=217)  2145    376        3 ( 597/ 258- 781/ 529)  79.7 E44-E35 S33-S11
     11  569  317 (S24E54 L=202)   517    122        2 ( 546/ 279- 593/ 359)  50.2 E55-E50 S27-S21
  

    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