Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  531 1523 (N32E33 L=326)   935    271        1 ( 478/1379- 594/1631)  73.6 E34-E29 N22-N39
      2 1564 1343 (N14W37 L= 36)   281     72        1 (1531/1318-1591/1368) -53.1 W34-W39 N12-N16
      3 1288 1277 (N12W17 L= 17)   230     47        1 (1272/1250-1308/1312) -90.0 W16-W18 N09-N13
      4  845 1157 (N07E10 L=349)   663    229        1 ( 746/1062- 906/1227) -54.0 E15-E06 N00-N11
      5  463 1099 (N06E35 L=324)   314     69        1 ( 433/1077- 499/1124)   0.0 E37-E32 N04-N07
      6  461 1037 (N02E35 L=324)   520    139        1 ( 426/1007- 508/1088) -38.7 E37-E32 N00-N05
      7 1787  685 (S26W55 L= 54)   680    160        0 (1749/ 621-1832/ 772)  63.4 W52-W58 S30-S20
      8 1408  508 (S37W25 L= 25)   697    151        2 (1330/ 484-1481/ 542)   0.0 W19-W31 S37-S34
      9 1225  390 (S44W12 L= 11)   597    151        0 (1142/ 369-1303/ 424)   8.7 W05-W18 S46-S42
  

    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