Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  569 1737 (N24E55 L=303)   635    183        0 ( 524/1708- 613/1763)  -4.8 E61-E49 N22-N25
      2  282 1189 (S15E53 L=305)   591    126        0 ( 266/1132- 297/1262) -63.4 E55-E50 S18-S10
      3 1162 1198 (N07W03 L=  1)   381    118        1 (1120/1154-1186/1262)  90.0 W01-W04 N03-N11
      4  748 1151 (S07E19 L=339)   326     85        1 ( 714/1111- 789/1176)  53.1 E20-E16 S08-S04
      5 1795 1085 (N21W47 L= 45)  1040    206        0 (1743/ 979-1841/1157) -32.5 W41-W54 N16-N23
      6  271 1056 (S23E51 L=307)   919    215        1 ( 205/1026- 348/1093)   0.0 E57-E43 S24-S21
      7  908  756 (S24W01 L=359)   342     81        0 ( 861/ 737- 947/ 778)  59.0 E01-W03 S26-S21
      8  883  627 (S32W04 L=  2)   680    183        1 ( 802/ 593- 973/ 685)   0.0 E03-W10 S34-S30
      9 1144  520 (S30W24 L= 22)   375     82        1 (1112/ 500-1186/ 556)  26.6 W21-W26 S32-S28
     10 1700  480 (S14W66 L= 63)   737     79        0 (1663/ 442-1735/ 525)  71.6 W63-W67 S16-S10
     11 1287  317 (S37W45 L= 43)   591    128        1 (1227/ 301-1356/ 333)  40.6 W41-W48 S39-S33
  

    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