Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  406 1429 (N30E42 L= 16)  3552    520        8 ( 267/1227- 524/1618) -32.0 E78-E30 N15-N45
      2 1581 1488 (N22W44 L=102)   762    185        2 (1491/1461-1658/1511)  -8.1 W37-W51 N20-N24
      3  794 1415 (N25E11 L= 47)   214     97        2 ( 758/1374- 826/1449) -56.3 E13-E09 N21-N27
      4  188 1235 (N20E62 L=356)  2240    148        6 ( 136/1183- 247/1294)  28.6 E70-E55 N16-N23
      5  128 1084 (N11E69 L=349)  2369    276        8 (  87/ 997- 190/1182) -42.7 E78-E60 N05-N17
      6  452  703 (S15E41 L= 16)   861    273        0 ( 377/ 605- 527/ 805)  34.5 E50-E33 S20-S09
      7  253  603 (S19E66 L=352)  2614    336        5 ( 214/ 468- 302/ 771)  59.9 E73-E58 S27-S08
      8 1293  743 (S20W14 L= 72)   275    103        0 (1255/ 711-1342/ 765)  26.6 W11-W17 S22-S19
      9 1109  675 (S23W02 L= 60)   779    279        1 ( 999/ 641-1208/ 751) -36.0 E04-W08 S25-S17
     10 1490  465 (S41W32 L= 90)   655    110        2 (1448/ 437-1530/ 504)  -7.1 W27-W35 S43-S38
  

    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