Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  498 1655 (N36E51 L=162)   680    117        0 ( 456/1636- 553/1672)  20.6 E54-E46 N34-N37
      2  603 1651 (N37E42 L=171)   444     97        1 ( 575/1615- 631/1671)  90.0 E45-E39 N34-N38
      3  770 1299 (N16E19 L=194)   288    116        0 ( 736/1258- 821/1319)  31.0 E21-E15 N12-N16
      4  209 1223 (N05E64 L=150)   586     66        0 ( 193/1186- 231/1260)  68.2 E65-E61 N02-N07
      5  656 1135 (N04E24 L=189)   308     59        1 ( 632/1103- 677/1167)  63.4 E25-E22 N02-N06
      6 1091 1110 (N07W03 L=217)   229     75        0 (1054/1082-1132/1138) -26.6 W00-W06 N05-N08
      7  409 1053 (S03E41 L=173)   245     64        1 ( 395/1025- 428/1087) -56.3 E42-E39 S04-S01
      8 1098  889 (S06W06 L=219)   749    221        1 (1056/ 769-1142/ 983) -61.7 W02-W09 S13-N00
      9  728  871 (S11E17 L=197)   425    146        2 ( 651/ 853- 792/ 893)  11.3 E22-E12 S12-S09
     10 1186  682 (S18W14 L=227)   275     53        0 (1157/ 662-1217/ 699)  18.4 W12-W16 S19-S17
  

    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