Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1152 1131 (N10W10 L= 88)   457    142        1 (1130/1075-1171/1203)  60.3 W07-W12 N06-N14
      2  398  939 (N15E43 L= 36)   928    231        1 ( 337/ 845- 446/1024) -63.4 E46-E39 N08-N20
      3 1498  733 (S22W22 L=101)  2617    761        2 (1225/ 507-1687/ 975)  20.6 W00-W40 S30-S12
      4  221  885 (N14E61 L= 17)  1159    174        5 ( 188/ 811- 267/ 931)  49.4 E65-E55 N09-N16
      5 1265  708 (S17W07 L= 85)  1007    474        2 (1197/ 541-1346/ 881)  50.2 E01-W15 S26-S08
      6  377  612 (S04E54 L= 24)   677     83        1 ( 356/ 578- 411/ 653) -78.7 E56-E51 S06-S01
      7 1559  574 (S35W25 L=104)   540    117        0 (1517/ 526-1595/ 618)  35.5 W21-W28 S37-S32
      8  575  305 (S28E56 L= 22)   521     80        1 ( 544/ 289- 614/ 331) -21.8 E59-E52 S30-S25
      9  850  249 (S39E38 L= 40)  1195    182        3 ( 793/ 181- 893/ 326)  30.1 E51-E32 S43-S32
  

    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