Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1164 1758 (N43E18 L=197)   683    188        0 (1097/1682-1213/1825) -22.4 E27-E10 N38-N47
      2  315 1448 (N00E61 L=154)  1080     79        2 ( 291/1404- 345/1486) -78.7 E65-E56 S02-N03
      3 1161 1386 (N18E03 L=213)   722    171        0 (1093/1320-1221/1435) -18.4 E07-W02 N15-N19
      4  352 1349 (S05E52 L=163)   550     58        0 ( 328/1311- 377/1383)  90.0 E54-E50 S07-S02
      5  204 1072 (S24E61 L=155)   501     50        2 ( 182/1042- 223/1104) -76.0 E63-E58 S25-S21
      6  750  918 (S20E13 L=203)  1088    431        1 ( 565/ 823- 923/ 976)  40.4 E23-E03 S29-S12
      7  387  812 (S36E38 L=177)   712    188        2 ( 357/ 728- 446/ 873)  74.7 E40-E32 S41-S30
      8  881  762 (S25E01 L=215)   562    181        0 ( 809/ 717- 968/ 793)  56.3 E03-W03 S29-S20
      9  960  635 (S30W08 L=224)   951    299        3 ( 815/ 614-1085/ 664)  28.1 W00-W15 S34-S25
     10 1260  505 (S28W32 L=247)   480    131        1 (1202/ 485-1318/ 532)  12.5 W27-W36 S28-S25
     11 1552  468 (S20W54 L=270)   951    172        2 (1482/ 425-1628/ 509)  63.4 W51-W58 S24-S14
     12  375  460 (S62E48 L=168)   976    115        1 ( 319/ 434- 427/ 492)   2.6 E60-E36 S64-S59
  

    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