Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1755 1427 (N42W40 L= 77)   420     97        1 (1730/1406-1777/1454)  45.0 W37-W43 N39-N44
      2 1363 1347 (N21W10 L= 47)   315    127        0 (1320/1284-1409/1400) -32.0 W06-W14 N18-N23
      3 1182 1327 (N15W00 L= 37)   217     59        1 (1150/1306-1208/1356)  76.0 E01-W01 N13-N17
      4  260 1285 (S10E59 L=338)   385     88        2 ( 244/1246- 288/1325)  90.0 E61-E56 S12-S07
      5 1762 1082 (N19W44 L= 81)  1038    192        1 (1739/ 988-1783/1182) -55.0 W40-W48 N14-N25
      6 1886  914 (N15W62 L= 99)  1799    204        5 (1860/ 842-1908/1003) -37.9 W56-W66 N11-N18
      7 1223  953 (S05W13 L= 49)   355     97        1 (1178/ 926-1253/ 995)  53.1 W10-W14 S06-S02
      8  548  813 (S31E24 L= 13)  1229    431        0 ( 369/ 762- 730/ 877)   8.1 E38-E10 S33-S27
      9 1380  759 (S10W28 L= 64)   416    113        0 (1343/ 704-1429/ 816)  83.7 W26-W29 S14-S05
     10 1264  604 (S22W27 L= 64)   479     85        3 (1206/ 578-1310/ 634)  45.0 W24-W29 S24-S19
     11 1625  347 (S23W68 L=104)   963    126        1 (1570/ 307-1679/ 394)  63.4 W65-W70 S26-S18
  

    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