Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1019 1636 (N39E12 L=124)   667    166        4 ( 989/1576-1064/1689) -53.1 E14-E07 N35-N43
      2 1419 1560 (N41W21 L=157)   563    143        1 (1373/1519-1465/1599)  41.2 W17-W25 N37-N44
      3  400 1451 (N17E51 L= 84)   527    128        1 ( 352/1420- 437/1498)  39.8 E55-E48 N14-N20
      4 1323 1463 (N32W14 L=149)   284     45        1 (1301/1437-1345/1486) -56.3 W12-W15 N30-N33
      5  760 1346 (N16E22 L=114)   228     83        1 ( 709/1328- 797/1364)  31.0 E24-E19 N13-N17
      6  966 1247 (N13E07 L=129)   225     89        0 ( 921/1232-1011/1262)  21.8 E09-E04 N11-N14
      7  989 1151 (N07E04 L=132)   763    305        3 ( 876/1053-1092/1233) -25.0 E12-W03 N02-N10
      8  631  945 (S10E23 L=112)   314    114        0 ( 578/ 926- 683/ 976)  33.7 E26-E20 S11-S07
      9  984  855 (S11W00 L=136)   711    155        1 ( 953/ 772-1009/ 926) -71.6 E02-W02 S15-S06
     10  440  567 (S38E38 L= 98)   391    115        0 ( 406/ 540- 480/ 600)  32.0 E41-E33 S40-S35
  

    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