Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  935 1391 (N11E14 L=243)  2000    503        3 ( 838/1216-1084/1608) -35.8 E26-E01 N04-N22
      2 1148 1544 (N26E06 L=251)   252     92        0 (1130/1495-1161/1586) -51.3 E08-E04 N23-N28
      3  346 1433 (N02E55 L=202)  1975    451        5 ( 229/1335- 437/1478)  14.7 E68-E46 S02-N05
      4 1495 1329 (N23W20 L=277)   502     87        1 (1470/1281-1515/1377) -71.6 W18-W22 N19-N25
      5  208 1111 (S19E58 L=199)  2288    473        6 ( 159/1030- 259/1218) -48.0 E64-E51 S22-S12
      6 1488 1067 (N07W26 L=283)  1205    238        2 (1449/ 947-1528/1187) -47.7 W20-W30 N01-N13
      7  652  867 (S25E19 L=238)   814    268        1 ( 531/ 832- 774/ 917)  31.0 E26-E11 S28-S19
      8  260  786 (S38E50 L=207)  1618    426        3 ( 170/ 728- 385/ 820)  22.2 E60-E38 S42-S33
      9 1640  410 (S22W63 L=320)  1748    181        2 (1571/ 368-1691/ 476)  66.0 W58-W66 S26-S17
  

    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