Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1311 1717 (N36W43 L=102)   578    104        0 (1273/1690-1353/1735) -18.4 W39-W46 N33-N37
      2 1611 1493 (N14W51 L=109)   400     80        0 (1577/1462-1642/1530) -68.2 W49-W53 N11-N16
      3  985 1425 (N28W09 L= 67)   226     62        0 ( 955/1401-1021/1440)   0.0 W06-W11 N27-N29
      4 1084 1263 (N16W10 L= 68)   235     97        0 (1046/1232-1116/1306)  31.0 W07-W12 N14-N18
      5  149 1253 (N39E61 L=358)   757    210        1 ( 116/1233- 185/1271) -20.6 E64-E56 N37-N40
      6  826 1201 (N20E07 L= 51)   401    117        1 ( 761/1181- 879/1221) -29.7 E10-E03 N17-N21
      7 1463 1004 (S08W24 L= 82)   206     37        1 (1442/ 982-1488/1029) -18.4 W22-W25 S09-S06
      8  267  963 (N19E51 L=  7)   443    127        0 ( 226/ 939- 311/ 980)  -7.1 E55-E47 N17-N19
      9 1798  734 (S35W46 L=104)  2043    188        4 (1726/ 632-1847/ 848)  45.0 W38-W51 S39-S28
     10 1303  601 (S27W05 L= 63)  1101    354        1 (1274/ 442-1349/ 735)  70.7 W00-W11 S38-S18
     11  962  508 (S23E18 L= 40)   855    295        0 ( 888/ 408-1007/ 643)  36.3 E26-E11 S27-S15
     12  679  465 (S18E37 L= 21)   836    195        3 ( 641/ 411- 725/ 539)  69.4 E39-E33 S22-S14
  

    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