Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1171 1857 (N65E08 L= 82)  5204    455        5 ( 978/1777-1361/1931)  -1.1 E44-W18 N58-N72
      2 1128 1681 (N45E05 L= 85)  1321    193        3 (1042/1626-1202/1724) -11.3 E13-W02 N41-N47
      3  531 1667 (N33E54 L= 36)  1006    237        3 ( 472/1590- 627/1707)  70.0 E58-E46 N26-N37
      4  342 1507 (N19E62 L= 28)   493     76        2 ( 323/1482- 373/1542) -20.6 E65-E56 N17-N21
      5  328 1358 (N10E55 L= 34)  1059    184        5 ( 289/1287- 382/1430) -36.0 E62-E49 N05-N13
      6  692 1385 (N17E27 L= 63)   247     82        0 ( 670/1352- 723/1422)  78.7 E28-E26 N14-N19
      7 1507 1196 (N18W29 L=119)  1235    346        1 (1400/1079-1593/1311) -36.9 W20-W36 N11-N23
      8  442 1150 (S01E40 L= 50)   240     65        0 ( 426/1132- 470/1179) -14.0 E41-E37 S02-N00
      9 1857 1077 (N16W62 L=152)  1145    160        0 (1839/1014-1882/1154) -69.4 W60-W66 N12-N20
     10  344  987 (S12E46 L= 44)   286    101        1 ( 308/ 950- 376/1030)  51.3 E48-E43 S14-S09
     11  153  862 (S23E71 L= 19)   943     73        3 ( 130/ 835- 172/ 895) -36.9 E77-E66 S24-S21
     12  261  791 (S26E55 L= 35)   382     83        0 ( 243/ 748- 279/ 831)  45.0 E58-E52 S29-S23
     13 1203  736 (S15W16 L=106)   248     60        0 (1173/ 713-1237/ 756)  45.0 W13-W17 S16-S13
     14  507  701 (S28E31 L= 59)  1371    316        2 ( 349/ 658- 615/ 741)  20.6 E46-E22 S33-S23
     15 1203  609 (S23W18 L=108)   334    136        1 (1165/ 553-1249/ 661) -33.7 W14-W23 S25-S19
     16 1040  595 (S26W08 L= 98)   382    154        0 ( 987/ 550-1099/ 651)  60.9 W05-W10 S30-S21
     17 1435  508 (S26W38 L=128)   296     44        1 (1419/ 482-1457/ 534) -63.4 W36-W40 S27-S23
     18  675  218 (S69E26 L= 64)   941    100        2 ( 648/ 181- 704/ 260)  71.6 E32-E18 S75-S62
  

    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