Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  789 1578 (N19E29 L=169)   409    148        0 ( 767/1514- 804/1651) -60.9 E32-E26 N15-N24
      2  986 1501 (N19E14 L=183)   389    168        0 ( 933/1430-1029/1565) -35.0 E19-E09 N15-N22
      3  566 1451 (N07E38 L=159)   321    110        0 ( 516/1430- 611/1471)  38.7 E41-E35 N05-N09
      4  882 1424 (N12E18 L=179)   244    101        0 ( 869/1371- 896/1471) -74.1 E19-E16 N08-N15
      5 1345 1426 (N24W09 L=206)   230     87        0 (1328/1384-1362/1463) -53.1 W07-W10 N21-N26
      6 1377  984 (S01W21 L=218)   370    148        0 (1302/ 957-1443/1005)   6.3 W16-W25 S02-N00
      7 1102  972 (S08W06 L=203)   181     90        0 (1066/ 957-1150/ 997)  11.3 W03-W08 S08-S07
      8  256  827 (S35E54 L=144)   564    177        0 ( 218/ 796- 296/ 872) -29.7 E58-E48 S37-S31
  

    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