Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  858 1650 (N50W05 L= 62)   305     69        1 ( 832/1623- 894/1674)  15.9 W02-W09 N48-N51
      2 1835 1166 (S05W60 L=118)   940    114        3 (1800/1115-1855/1214)  80.5 W57-W63 S08-S02
      3  727 1109 (N18E16 L= 41)   267     72        1 ( 696/1085- 760/1137)   0.0 E18-E13 N16-N19
      4  266  903 (N13E55 L=  2)  1294    162        4 ( 237/ 836- 315/ 987) -77.5 E58-E51 N08-N17
      5 1255  453 (S34W02 L= 60)  1533    357        4 (1167/ 331-1342/ 598)  45.0 E07-W10 S41-S25
      6  938  447 (S27E19 L= 38)   555    164        1 ( 885/ 393- 970/ 520)  37.9 E24-E15 S29-S22
      7  557  417 (S22E48 L= 10)   591    108        2 ( 526/ 375- 581/ 477)  37.9 E52-E42 S24-S17
      8 1060  240 (S47E21 L= 37)  1185    181        4 ( 990/ 185-1121/ 300)  11.3 E29-E13 S51-S41
  

    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